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ABSTRACT 

In order to study how outflows from protostars influence the physical and 
chemical conditions of the parent molecular cloud, we have observed Barnard 1 
(Bl) main core, which harbors four Class and three Class I sources, in the CO 
(J = 1 - 0), CH3OH {Jk = 2k- Ik), and the SiO (J = 1 - 0) lines using the 
Nobeyama 45 m telescope. We have identified three CO outflows in this region; 
one is an elongated (~ 0.3 pc) bipolar outflow from a Class protostar Bl-c in 
the submillimeter clump SMM 2, another is a rather compact (~ 0.1 pc) outflow 
from a Class I protostar Bl IRS in the clump SMM 6, and the other is an extended 
outflow from a Class I protostar in SMM 11. In the western lobe of the SMM 
2 outflow, both the SiO and CH3OH lines show broad redshifted wings with the 
terminal velocities of 25 km s~^ and 13 km s~^, respectively. It is likely that the 
shocks caused by the interaction between the outflow and ambient gas enhance 
the abundance of SiO and CH3OH in the gas phase. The total energy input rate 
by the outflows (1.1 x IO^'^Lq) is smaller than the energy loss rate (8.5 x 10^'^Lq) 
through the turbulence decay in Bl main core, which suggests that the outflows 
can not sustain the turbulence in this region. Since the outflows are energetic 
enough to compensate the dissipating turbulence energy in the neighboring, more 
evolved star forming region NGC 1333, we suggest that the turbulence energy 
balance depends on the evolutionary state of the star formation in molecular 
clouds. 

Subject headings: stars: formation — stars: low mass — stars: pre-main sequence 
— ISM: clouds — ISM: jets and outflows — ISM: individual (Barnard 1) 



1. INTRODUCTION 

Molecular outflows are common features accompanied with young stellar objects (YSOs; 



Bachillerlll996[ ). The outflows have large impacts on the surrounding molecular cloud both in 
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the chemical and physical aspect. Shocks excited by the propagation of the outflow raise the 
temperature of surroundir ig material and liberate various molecules such as NH.^, CH.sOH 



and SiO from dust gra ins (lAUamandola et al.lll992l : ISchilke et al.lll997l : iBachiller et al.ll2001 



J0rgensen et al.ll2004l ). As a result, the relative abundance of these molecules against H2 
reaches several orders of magnitude larger than that in quiescent molecular clouds. 

The supersonic turbulence in molecular clouds is ubi quitous but should qu ickly dissipate 
through shocks, thus some driving sources are needed (IMac Low et al.lll998l ). Protostellar 
outflows are consi dered to be possible sources to maintain the supersonic turbulence in 
molecular clouds ( iNorman fc Silki Il980l ) and regulate the star formation, since the turbu- 
lence serves as a counterwork against the gravitati onal collapse in molecula r clouds. The 
observ ational studies of Herbig-Haro (HH) objects (jWalawender et al.ll2005al ) and CO out- 
flows (IKnee &: Sandelll I2OOOI : IStanke fc Williams! 120071 ) compared the energy input by the 
outflows and turbulence energy decay, and found that the outflows have enough energy to 
retain the turbulence in the molecular clouds in at least local scale (~ 0.5 pc), such as NGC 
1333 or L1641-N, but not whole Perseus or Orion giant molecular cloud scale (> 10 pc). The 
theoretical simulation also pointed out that the protostellar outflow s play an important r ole 
to retain the turbulence and regulate the star formation efficiency ( iLi fc Nakamural 120061 ). 



Barnard 1 (Bl) is one of the star-forming regions locat ed in the Perseus raolecul ar cloud 
complex. The distance to Bl is estimated to be 230 pc by ICernis fc Straizysl ( 120031 ) on the 
basis of the extinction study. This value is consistent with the distance of 235 ± 15 pc to the 
NGC 1333 molecular cloud, which is located at ~ 1° west of Bl, derived from the p arallax 
measurements of the H2O masers associated with the YSO SVS 13 (IHirota et al.ll2008l ). Note 
that this distance is different from the previously adopted values of 350 pc and 318 pc, which 
are th e distances to the Perseus 0B2 association ( Borgman fc Blaauw 1964 de Zeeuw et al. 



19991 ). Bl was observed in several molecular lines by IBachiller et al.l ( 1l990l ). They revealed 
that the dense part of the cloud traced by the C^^O (1 — 0) line has a size of ~ 1.3 x 3.5 
pc and a mass of 518 Mq. Inside this region, a denser portion called "main core" with a 
diameter of 0.5 pc and a mass of 65 Mfri □ w as found in the NH3 emission. A submillimeter 
continuum survey ( iWalawender et al.ll2005bl ) identifled 12 submillimeter clumps in Bl and 
six of them are located in the main core region. Two of the submillimeter clumps (SMM 1 
and SMM 2) harbor Class YSOs and one (SMM 6) harbor s a Class I YSO. I nside SMM 
1, there are two extremely young sources, Bl-bN and Bl-bS (IHirano et al.lll999l ). A Class I 
source in SMM 6 was detect ed by IRAS and is re f erred as IRAS 03301-1-3057 or Bl IRS, and 
accompanied by an outflow (IBachiller et al.lll990l : lYamamoto et al.lll992l : iHirano et al.lll997l : 
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de Gregorio-Monsalvo et al.ll2005l ). The YSO in SMM 2 (also referred as Bl-c) located at 



~ 1' north of IRAS 03301+3057 is also a driving source of the outflow that is extending along 
the east-west direction. In addition to those four YSOs (Bl IRS, Bl-bN, Bl- bS, and Bl-c), 



one Cl ass YSO was identified in SMM 3, and one Class I source in SMM 11 (IHatchell et al. 



2007al ). There is another Class I source, LkHa 327, without any submillimeter counterpart. 



In total, four Class and three Class I sources are identified in Bl main core. The properties 
of the YSOs in Bl main core are summarized in Tabled! Near infrared observations in Ha 
and H2 rovibrational lines found a number of shock s produ ced by protostellar outflows in 

2005al jbl. 



the Bl main core region (iWalawender et al 



20091 ). The most prominent shocks 



are associated with the outflow from SMM 2. Mid-infrared observations with Spitzer Space 
Telescope have clearly revealed the S-shaped features of the knotty chains of the shocks 
cen tered at SMM 2, as w ell as several shocked regions around SMM 1, SMM 6, and SMM 
11 ( |j0rgensen et al.ll2006[ ). 

This paper presents multi-molecular line observations of Bl main core to investigate the 
effects of protostellar outflows on the physical and chemical properties of the surrounding 
core. The details of our observations are described in Section 2. The results of the observa- 
tions in CO, CH3OH, and SiO hues are described in Section 3, followed by the discussion on 
the chemical abundance of the molecules and turbulence energy budget in Section 4 and a 
summary in Section 5. 



2. OBSERVATIONS 

CO (J = 1 - 0), CH3OH (Jk = 2k - Ik), and SiO (J = 1 - 0) line observations 
were performed in 2006 April with the 45 m telescope at Nobeyama Radio Observatory 
(NRO). The frequencies of the observed lines a re summarized in Table [21 The 25-Beam Array 



Receiver System (BEARS; ISunada et al.ll2000l ) operated in the double sideband (DSB) mode 



was used for the CO and CH3OH observations. The system noise temperature of BEARS 
during the observations was in the range of Tgys (DSB) = 200 - 500 K. The main beam 
efficiency was ?7mb ~ 0.45 at 115 GHz. The half power beamwidth (HPBW) was 18"at 97 
GHz (CH3OH) and 15"at 115 GHz (CO), and the separation between each beam was 41'.'1. 
Twelve pointings of BEARS were adopted to map the 10' x 6^5 region, which covers the 
entire region of Bl main core, with a grid spacing of 29'.'1. The intensity scales of the 25 
receiver beams of BEARS at 97 GHz were calibrated by comparing the intensity of NGC 
2264 observed by each beam of BEARS with that measured by the single-beam receiver SIOO. 
The intensity scales of the 25 beams at 115 GHz were calibrated using the scaling factor table 
provided by the observatory. The back end was an autocorrelator with a bandwidth of 32 
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MHz and a frequency resolution of 37.8 kHz, which corresponds to a velocity resolution of 
0.098 km s"^ and 0.12 km s'^ at 115 GHz (CO) and 97 GHz (CH3OH), respectively. 

The SiO {J = 1 — 0) line was observed along the jets from SMM 2 and SMM 6 with 
a single-beam SIS receiver S40 operated in the single-sideband (SSB) mode. The position 
angles of the strip scans are 77° and 133° for SMM 2 and SMM 6, respectively. The obser- 
vational grid spacing was 38'.'4, which was the same as the HPBW at 43 GHz. The system 
noise temperature Tgyg (SSB) during the observations was ~ 180 — 300 K. The main beam 
efficiency was t^mb = 0.77 at 43 GHz. The bandwidth of the autocorrelator was set to 32 
MHz. The spectral resolution was 37.8 kHz which corresponds to a velocity resolution of 
0.26 km s-i at 43 GHz. 

The telescope pointing was checked by observing the SiO maser emission from NML 
Tau every 60-90 minutes and the error was estimated to be less tha n ^ 6" in 100 GHz and 



10" in 43 GHz observations. The standard chopper-wheel method (IKutner fc Ulichlll98ll ) 
was adopted to convert the received signal into the antenna temperature. The data were 
reduced using NEWSTAR, a data reduction software developed by NRO. Linear baselines 
were subtracted in most cases. 



3. RESULTS 
3.1. Molecular Outflow Traced by the CO Line 

3.1.1. Distribution of the Outflows 

Figure [T] displays the CO (1 — 0) line profiles observed at four representative positions 
in Bl main core, i.e. (a) 120" west of SMM 2, (b) SMM 6, (c) SMM 1, and (d) 5f5 southwest 
of SMM 6. In this paper, these four positions are referred to as positions "a", "b", "c", and 
"d" , re spectively. Since there is no jet-like component around this position in the infrared 



image (jWalawender et al.ll2005bl ). it is likely that the spectrum observed at the position d 
is not affected by outflows. Therefore, the velocity range of the ambient cloud component, 
"line core", was determined to be from Vlsr = 2.1 km s~^ to Vi,sk = 9.6 km s~^, using the 
CO spectrum at the position "d" . The CO emission outside of this velocity range is regarded 
as the outflow component. The CO spectra obtained at the positions "a" and "b" exhibit 
line wings in the redshifted side and blueshifted side, respectively. The most prominent CO 
redshifted wing was observed at the position "a" . At the position "c" , wing-like blueshifted 
emission is detected in the velocity range of 1 km s~^ < Vlsr < 3 km s~^. This component 
might be originated from an outflow associated with SMM 1, however, because of its low 
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velocity, it is unclear whether this component arises from the outflow or not. In addition, 
the spatial distribution of the blueshifted CO wing (see below) is not localized on SMM 1. 
Therefore, in this paper we do not include this component as an outflow candidate. 

The spatial distributions of the blue- and redshifted outflow components are shown in 
Figj2l Although the spatial distributions of the high velocity components are complex, two 
bipolar outflows centered at SMM 2 and SMM 6, and the extended blueshifted emission at the 
southeast corner are recognized. The outflow centered at SMM 2 extends along the east— west 
direction. The extent of the westw ard redshifted component is 120"(0.13 pc) and reaches the 



position of a shocked region MH 2 ( iWalawender et al 



2005bl ) . The redshifted emission peaks 



at the head of the lobe (position "a" ) and the velocity is also the largest at this position. On 
the other hand, the eastern blueshifted component peaks at 150"(0.17 pc) southeast of SMM 
2. Neither the redshifted peak in the west nor the blue shifted peak in the east was covered 
by th e previous CO (3 — 2) map (IHatchell et al.ll2007bl ) and the interferometric CO (1 — 0) 
map ( iMatthews et al.l 120061 ). Hence the entire picture of this outflow has been revealed by 
our observations for the flrst time. There is a compact blueshifted component in the western 
lobe and an extended redshifted component in the eastern lobe. This suggests that the 
outflow axis is close to the plane of the sky, and the overlapping blueshifted and redshifted 
components are due to the projection of the front and back sides of the lobes. The outflow 
traced by the CO line is well ali gned with the S-shaped knots seen in the Spitzer IRAC band 
2 image ( jJorgensen et al.l 120061 ). This IRAC band contains several H2 pure-rotational lines 
that are considered to arise from shocked regions. Therefore, it is likely that the S-shaped 
feature traces the shocks in the outflow driven by SMM 2. 

The outflow from SMM 6 is more compact than that from SMM 2. The blueshifted 
emission is detected toward the position of SMM 6 and slightly extended (~ 40", 0.04 pc) 
to the northeast. There is a redshifted component which peaks at 65" (0.07 pc) southwest 
of SMM 6. It is uncertain whether this redshifted component is originated from the SMM 
6 outflow or not, because there is another YSO candidate SSTc2d J033314. 4+310711 at 20" 
(0.02 pc) east of the redshifted peak. This source is classifled as a "YSOc_red" (YSO candi- 
date and red) in the Spitzer c2d catalog, and might be the driving source of the redshifted 
outflow. However, the grid spacing of the current observations is not flne enough to pinpoint 
the position of the driving source unambiguously. Therefore, in the following discussion, the 
redshifted component at 65" southwest of SMM 6 is treated as a redshifted lobe of the SMM 
6 outflow. A series of shocke d knots MH 4 was found at 15" west of SMM 6 and thought 
to be emanated from SMM 6 (IWalawender et al.ll2005bl . l2009l ). These knots align east-west 
direction and are not parallel to the SMM 6 CO outflow. Since the length of the knots is 
20", smaller than our mapping grid spacing, the relation between the knots and the SMM 6 
outflow is not clear in the present observations. An elongated reflection nebula was identi- 
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fied in the i^s-band image around SMM 6 (IWalawender et al.ll2005bl ) and the orientation is 
consistent with that of the SMM 6 outflow. This consistency would support the idea that 
SMM 6 is the driving source of the redshifted component southwest of SMM 6. 

The extended blueshifted CO emission is found around the shocked knots MH 8 and 
HH 790. In the CO image, this blueshifted component is continuously connected to the 
blueshifted lobe of the SMM 2 outflow. If this blueshifted component in the southeast area 
is a part of the SMM 2 outflow, however, the outflow is highly asymmetric with a blueshifted 
lobe more than twice as longer as the redshifted lobe. In addition, the blueshifted lobe 
bends sharply to the south at the position around the MH 1 knot. In contrast, the shocked 
components seen in the Spitzer image suggest that the eastern lobe and the western lobe 
of the SMM 2 outflow are similar in length. Therefore, it is possible that the blueshifted 
emission in the southeastern area originates from the other outflow. One possible candidate 
is the giant east-west outflow from SMM 11, which was identified in the deep H2 image 
obtained by IWalawender et al.l (120091 ). This outflow is parallel to the SMM 2 outflow, and 
contains MH 8 and HH 790 in the eastern lobe. If MH 8 and HH 790 are related to the 
outflow from SMM 11, it is natural to consider that the blueshifted CO component around 
MH 8 and HH 790 is also a part of the SMM 11 outflow. In the H2 image, the western 
lobe of the SMM 11 outflow extends to MH 10 and MH 3. The redshifted CO emission was 
detected around MH 10 and would be the part of the SMM 11 outflow. However, the giant 
outflow extending further west was not clearly detected in our CO observations. 

The outflows from SMM 1, SMM 3, and LkHa 3 2 7 are ii ot clearly seen in our CO (1 — 0) 
map, although they were reported by iHatchell et al.l (l2007bl ) on the basis of the CO (3 — 2) 
observations, except for that from LkHa 327. Our CO (1 — 0) maps show the redshifted 
er nission to the northeas t of SMM 1. This component is also seen in the CO (3 — 2) map 
of IHatchell et al.l ( l2007bl ). and identified as the redshifted lobe of the SMM 1 outflow. It is 
possible that this red shifted emission is related to the outflow from SMM 1 as argued by 
Hatchell et al.l ( l2007bl ). However, we regard this emission as a part of the outflow from SMM 

11 because there is no clear blueshifted component and bipolarity around SMM 1. The mid- 
and near infrared images show a chain of shocked knots that extends from SMM 1 to SMM 
3. The axis of this knot complex is parallel to that of the SMM 6 outflow. Although SM M 1 
is designated to be a driving source of these shocked features (IWalawender et al.ll2009[ ). no 
corresponding CO high-velocity component was detected in the present observations. LkHa 
327 was classified as Class I based on the s pectral energy distribution slope in the IRAS 

12 — 100 fim data (IWalawender et al.ll2005bl ). On the other hand, this source was classified 
as Class II based on the Two Micron All Sky Survey if-band, Spitzer IRAC, and MIPS data. 
The facts that LkHa 327 is located at the periphery of Bl main core and not associated 
with any 850/im emission indicate that the source is in the rather evolved stage. 
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3.1.2. Estimation of the Kinematic Properties of the Outflows 
The mass, momentum, kinetic energy, and th e momentum flux o f the identifled CO 



outflows are calculated with the method described by lBourke et al.l (Il997l ). In this procedure, 
we divided the high-velocity CO emissions originated from the outflows into the 0.4 km s~^ 
velocity bins, and calculated the quantities on an assumption of the local thermal equilibrium 
(LTE) condition. First, the column density of the molecule in the velocity range of v to 
V + Af , N{v), is given as 

ShZjT) exp{hBJiJ+l)/kT,^} ^ r+^" 
"^^"^ = 8vr3;.^(J+l) l-exp(-WA:T.- ^1 



'-ex J J V 



where 

f = \ (2\ 

^ J(Tex)- J(Tbg)l-exp(-r)' 
Tirp\ hvjk 
^^^^ = exp(/iz/AT)-l' 

h is the Planck constant, Z{T^ is the partition function, [i is the dipole moment of the 
molecule, B is the rotation constant of the molecule, k is the Boltzmann constant, v is 
the frequency of the line, r is the optical depth, Tgx is the excitation temperature, and Tbg 
is the background temperature. For linear molecules such as CO, the partition function 
is approximated to be 2^(T) = kT/hB if hB ^ kT. In the case of the CO (1 — 0) line 
under the optically thin assumption, the above equations become the following as given by 



Bourke et all (119971 ): 



N{v) = 2.31 X 10^^- , , ' wrr / UT ^ / 

1 - exp(-5.53/Tex) J(Tcx) - J{Tbg) Jv 

The mass at each observed position (a, 6) and velocity v, M(t>, a, 6) is given as 

M{v, a, S) = [H2/C0]/i^mHA(a, 5)N{v), (5) 

where [H2/CO] is the reciprocal of the CO abundance and adopted to be 10^, /im = 2.72 is 
the mean molecular mass per H2, mn is the mass of a hydrogen atom, and A{a, 6) is the area 
of the CO emission at a given observed position. Since the observing grid spacing (29'.'1) is 
larger than the telescope beam (15"), we adopted the squared value of the observing grid 
spacing as A{a, 6) under the assumption of the uniform emission distribution in the 29'.'1^ 
region. The mass at each velocity bin M{v) is estimated by integrating M{v, a, 6) over the 
outflow spatial extension, and then the total mass is derived by integrating M{v) over the 
outflow velocity range. 
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Since the optical depth of the CO hne is not obtained in our observations, the assumption 
of the optically thin CO (1 — 0) line is adopted, hence the estimation should be a lower limit. 
The exci tation temperature T p^ is estimated by comparing the CO (1 — 0) spectra and (3 — 2) 



spectra (IHatchell et al.ll2007bl ) as follows. Under the LTE condition, the total column density 



N and the column density at a given rotational energy level J, Nj, is given as 

N = Nj^^ exp{hBJ{J + l)/kT,,}. (6) 

Nj is also expressed as 

87ri/3 2J+1 1 /-^ , 

c-^ (2J + 3)Aj+ij 1 - exp{-hiy/kTc:,) J 

in the optically thin condition, where Aj+i j is the Einstein A coefficient for the transition 
from the level J + 1 to J. The relation between the excitation temperature Tgx and the line 
intensity ratio of the CO (3 — 2) and (1 — 0) transition is derived from Equations ([6]) and 

lTus,^-2dv _ 3,,p(_i6.6/y^^ ^l -exp(-16.6/T..) 



/ TMB,i-odv ' 1 - exp(-5.53/Tcx) 

The integrated intensities of the high-velocity CO (3 — 2) line toward the S MM 2 westward 



redshifted component and SMM 6 were calculated from the spectra shown by lHatchell et al. 



( 1200 7bl ). The integrated intensities, line ratios and the derived excitation temperatures are 
summarized in Table [31 Although the observed positions in the two transitions are slightly 
different, the offset (< 5") is smaller than 1/3 of the beam size and we ignore the positional 
difference. The derived excitation temperatures toward the two positions, 46 K and 42 K, are 
close to each other. We note that the derivation of the excitation temperature is sensitive to 
the accuracy of the relative calibration between the two independent observations with the 
different instruments. For example, a 20% uncertainty of the ratio yields ± 20 K uncertainty 
of the excitation temperature. Here we adopt 50 K as a value of the excitation temperature 
for the calculations o f the outflow parameters. This value is the same as that used by 



Hatchell et al.l ( l2007bl ) to calculate the outflow properties in the Perseus molecular cloud. 



The momentum P, kinetic energy E]^, dynamical timescale td, and the outflow momen- 
tum flux Foutflow are described as 

M{v){v -vo)dv, (9) 

= M{v){v-vofdv, (10) 

U = R/v,^ = R/iP/M), (11) 

^outflow = P/td, (12) 
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where vq is the systemic velocity of 6.7, 6.5, and 6.8 km s~^ for the SMM 2, SMM 6, 
and SMM 11 outflows, respectively, measured by fitting a Gaussian function to the CH3OH 
spectra toward the sources (see Section 3.2). R is the length of the outflow and Vch = P/M is 
the characteristic velocity. The mass-loss rate is estimated as M = -Foutfiow/^^windj where Vwmd 
is the velocity of the primary wind ejected from the close vicinity of the protostar. Since 
we do not have any inf ormation on the wind velocity from our observations, we assume 
^wind = 100 km s~^ (e.g., iGiovanardi et al.lll992l ). 



The inclination angle of the outflow axis is also unclear, therefore we employed 57?3 
from the line of sight, which is the mean inclination angle on the assumption of the random 
outflow orientation. This assumption should be reasonable for the SMM 2 outflow because 
the interferometric observations revealed that the opening angle of the outflow cavity in 
the 20"scale is ~ 55° and the red and blues hifted lobes are not overlapped at the position 
of the driving source ( iMatthews et al.l 120061 ). which indicates that the outflow axis inclines 
at the larger angle than the opening angle. This interpretation does not conflict with the 
fact that the blue- and redshifted components are overlapped in the larger scale as shown 
in Figj2l because these overlapped bl ue- and redsh i fted c omponents are interpreted as the 
cavity wall. For the SMM 6 outflow, iHirano et al.l (119971) suggested that the config u ration 
of the SMM 6 outflow is in a pole-on geometry, while Ide Gregorio-Monsalvo et al.l (|2005[ ) 
suggested that the outflow from SMM 6 is inclined toward southwest. Since it is difficult to 
estimate the inclination angle from our coarse grid observations and there is no consensus on 
the inclination angle of the SMM 6 outflow, it is reasonable to adopt the mean inclination 
angle. Since there is no previous study about the inclination angle for the SMM 11 outflow, 
we also adopt the mean inclination angle of 57?3. 

The blueshifted lobes of the outflows from SMM 2 and SMM 11 are not clearly separated 
in our observations. Therefore, we assumed that the transverse extent of the SMM 2 outflow 
lobe is symmetric with respect to the shocked knots, and regarded the blueshifted component 
at decl. > 31°08' as the SMM 2 outflow and that at decl. < 31°08' as the SMM 11 
outflow. The redshifted lobe of the SMM 11 outflow is also confused with the eastern 
redshifted component of the SMM 2 outflow. We assumed that the northeastern half of the 
redshifted component at the east of SMM2 belongs to the SMM 2 outflow, because this part 
is overlapped with the blueshifted lobe of the SMM 2 outflow. The southwestern half of this 
redshifted component is regarded as the redshifted lobe of the SMM 11 outflow. The outflow 
parameters derived on the above assumptions are summarized in Table HI 

The mass, momentum, and kinetic energy of the SMM 2 outflow are almost comparable 
to those of the SMM 11 outflow and a factor of 5-10 larger than those of the SMM 6 outflow 
The properties of outflows from Class and Class I sources were compiled by lBontemps et al. 
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( 119961 ) ■ and it was found that the outflow momentum flux is 1 order of magnitude larger 
for Class sources. The mean momentum flux is < Fqq >= 5.7 x 10~^ M(7i km s~^ yr~ ^ 



and 3 .8 x 10 ^Mq km s ^ yr ^ for Class and Class I sources, respectively ( iBontemps et al. 



19961 ). Considering that they adopted an opacity correction of a factor of 3.5, the momentum 
fluxes of the outflows emanated from SMM 2 and SMM 6 are in reasonable agreement with 
the mean values for Class and I sources, respectively. On the other hand, the outflow from 
SMM 11 has approximately 1 order of magnitude larger momentum flux compared to the 
mean value of the outflows from the Class I sources, probably because the spatial extent of 
this outflow is large. 



The momentum fluxes for the outflows were estimated by iHatchell et al.l (l2007b( ) based 
on their CO (3 — 2) observations. Converting the assumed distance to Bl from 320 pc to 
230 pc, the momentum fluxes are Fco = 1-2 x 10"^ 9.1 x 10"^ and 3.2 x lO^^M© km s"^ 
j^-i ^Qj^ ^Yie SMM 2, SMM 6, and SMM 11 outflow, respectively. The momentum flux they 
derived for all the outflows is smaller than our result, and the difference reaches a factor of 
70 for the SMM 2 outflow. There are several reasons for these discrepancies. First, they 
did not take the correction of the inclination i int o account, and this co rrection makes the 
momentum flux larger by a factor of sin i / cos^ i (IBontemps et al.l Il996l ) . Our assumption 
of the inclination of 57?3 yields the factor of 2.9. Second, we regard both the westward 



and ea stward redshifted components as the outflow from SMM 2, however, iHatchell et al. 



( j2007bl ) assigned the eastward redshifted component to the SMM 1 (Bl-bS) outflow. This 
different identification causes another factor of 2 discrepancy in the estirn ation of the outflow 
momentum flux. Finally, the mapped region by IHatchell et al.l (l2007bl ) does not cover the 
entire extent of the outflow from SMM 2. In general, the outflow momentum flux is a robust 
parameter against incomplete mappings. If a considerable fraction of the total momentum is 
carried by components outside the mapping area, however, incomplete mapping observations 
might result in underestimation of the momentum flux. In fact, the promi nent redshifted 



compon ent around the position "a" is located outside of the mapped area by IHatchell et al. 
(l2007bh . 



3.2. CH3OH 

The CH3OH line profiles at three representative positions, i.e. (a) 120" west of SMM 2, 
(b) SMM 6, and (c) SMM 1, are displayed in Fig. [3J The frequency separations between the 
2o — lo A'^ and the 2_i — l_i E lines and the 2o — lo -E and the 2o — lo A'^ lines correspond 
to 9.8 km s~^ and 6.3 km s~^, respectively. The widest line profile is seen at the position 
"a" (Figl3]a) where the peak of the CO redshifted outflow is located. At this position, the 
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2o — lo and 2_i — l^i E lines are blended because of the broad redshifted wings. The 
terminal velocity of the 2_i — l_iE line reaches 13 km s~^ away from the systemic velocity. 
The position "a" is located close (~ 10") to the branching point of the shocked feature 
seen in the Spitzer image (see FigJ2]). Similar blended CH3OH profiles are seen in three 
observed points around the position "a" . Toward SMM 6, both the 2o — lo A'^ and 2_i — l_i 
E lines show a blueshifted wing (Figl3]6) which is also seen in the CO line (Fig{T]6). No 
wing components are detected toward SMM l(Figl3]c) where no CO outflow was identified. 
From these results we divide the CH3OH emission into three velocity ranges; the blueshifted 
component (Vlsr < 5.5 km s^^), the line core component (5.5 km s^^ < Vlsr < 7.5 km s~^), 
and the redshifted component (Vlsr > 7.5 km s~^). 

The integrated intensity maps of the CH3OH lines are shown in FigJH Since the dis- 
tribution of the 2q — 1q line is similar to that of the 2„i — l^i E line, and these two 
lines are blended around the position "a", we integrated all the 2o — lo A'^ and 2_i — l_i E 
emission over the velocity of —3.6 km s^^ < Vlsr < 20.4 km s~^ with the reference frequency 
of 96.73939 GHz, the rest frequency of 2_i — l_i E line, and made a combined map which 
is shown in FigJH Hereafter we call this map "A + E map". The CH3OH {Jk = 2o — lo 
_i — l_i E) lines were detected around the dusty clump complex which contains SMM 
1, SMM 2, SMM 3, and SMM 6. There are two peaks seen in the A + E map; one peak is 
located at the position "a" and the other peak is located around SMM 6. No obvious peaks 
are detected in the eastern side of the SMM 2 outflow and around the other protostars. 
The obse rved CH3OH distribution is different from those of the dust continuum a nd NH3 
emission ( iMatthews Sz Wilson! I2OO2I : IWalawender et al.l l2005bl : ISachiller et al.l Il990l ) . Since 
one CH3OH peak located at the position "a" coincides with the peak of the redshifted CO 
emission, we examined the contribution of the high- velocity component to the A + E map by 
comparing the spatial distributions of the blue- and redshifted CH3OH emission with that 
of the line core emission (Fig. The high- velocity components peak at the position "a" 
and SMM 6, and it is obvious that the A + E map is largely affected by the high- velocity 
components. The line core emission shows a single peak at SMM 6, and its distribution is 
different from the distribution of the NH3 emission which peaks around SMM 1. The dis- 
tribution of the CH3OH line core emission is also different from the distri bution of the du st 



continuum emission which peaks at the positions of the SMM sources ( IWalawender et al. 



2005bl ). The coincidence of the peaks of the line core component and the high velocity 



(blueshifted) component at SMM 6 indicates that the peak around SMM 6 seen in the A + E 
map arises from both the quiescent gas and the outflow-related component. The CH3OH 
abundance enhancement due to the outflow - core interaction is a possible reason for the 
intensity enhancement around the outflows. The abundance enhancement is discussed in 
Section 4.1. 
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The CH3OH 2o — lo -E transition is much weaker than the other two transitions and 
detected above the Sa noise level at eight out of the 300 observed positions. An example of 
the spectrum is shown in FigI3]c, which is detected at {a, (5)(J2000) = (3^33"^20!8, 31°07'25'/1), 
10" southwest of the Class protostar SMM 1 (Bl-bS). The integrated intensity map of 
the 2o — lo -E" transition is shown in FigJH with red contours. This line shows its intensity 
peak around SMM 1 and extension to SMM 2 and 6. The CH3OH 2o — Iq E emission 
peak at SMM 1 is not seen in the other two transitions, whereas the high density tracers 
such as the NH3 (1,1), H^'^CO^ (J = 1 — 0) liries, and 850 /xm continuum emission show the 
emiss ion peaks around SMM 1 (IBachiller et al.lll990l : iHirano et al.lll999l : iMatthews &: Wilson 



2OO2I ). The upper state energy of the 2q — 1q E transition is 12.2 K, which is slightly higher 
than those of the other two transitions, 4.6 and 7.0 K f or the 2_i — l_i E and 2q — Iq A'^ 
transition, respectively (obtained from CDMS catalog, iMuUer et al.ll200l(). Since CH3OH 
emission tends to be sub-thermally excited in star-forming regions ( Isuckle fc Fuller 2002), 
the difference of the distributions between the 2o — Iq E transition and the other lower- 
excitation transitions is likely to reflect the density structure of Bl main core, as the 2o — lo 
E line and the high-density tracers show the common peak at SMM 1. 



3.3. SiO Emission along the Outflows 

The SiO {J = 1 — 0) emission was detected at four positions out of 13 observed positions 
along the SMM 2 outflow (Figj6]). The line profiles at the representative positions are shown 
in FigJTl In the northwestern lobe of the SMM 2 outflow, the SiO emission is redshifted and 
its line width increases as the distance from SMM 2 increases (from SMM 2 Wl to SMM 2 
W3, see FigJ2])- The terminal velocity of the SiO emission reaches up to 25 km s^^ at the 
position of SMM 2 W3, which corresponds to the peak of the redshifted CO emission. In 
the southeastern lobe, the SiO emission was detected only at one position labeled SMM 2E. 
This position corresponds to the peak of the blueshifted CO emission. However, the SiO 
emission peaks at Vlsr = 7.7 km s~^, which is redshifted with respect to the cloud systemic 
velocity. At the position of SMM 2 itself, there was no detectable SiO emission. 

The SiO emission originated from the SMM 6 outflow was searched for at four positions 
along the outflow, and the emission was detected only at the position of SMM 6. As shown 
in FigJTJ the SiO line at SMM 6 shows the narrow component (AV^ = 1.3 km s~^) with 
its peak velocity at the systemic velocity of Vlsr = 6.4 km s~^, and a possible blueshifted 
component with t he te rminal velocity of ~ 5 km s~^ with respect to the systemic velocity. 



Yamamoto et al.l ( 119921 ) also reported the detection of the SiO J = 1 — and 2 — 1 emission 



lines at 7" northwest from the position of SMM 6. They pointed out that the line profile 
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consists of a narrow [Av = 1.6 km s~^) component and br oad (Av = 6.4 krn . s~^) c omponent. 
The redshifted wing of the broad component reported in lYamamoto et aLl (jl992[ ) is weaker 
than the blueshifted wing, hence the redshifted wing was not detected probably due to the 
low signal-to-noise ratio in our observations. 



4. DISCUSSION 
4.1. CH3OH and SiO Abundance 

4.1.1. CHsOH 

The column density values of the CH3OH molecule at the positions of the protostars were 
calculated by assuming the LTE condition using Equation (1). Since CH3OH is not a linear 



molecule, the partition function was calculated from the energy level diagram ( iXu fc Lovas 



19971 ). The excitation temperature was ass umed to be 12 K, w hich is the same value as 



the kinetic temperature derived from NH3 ( iBachiller et al.lll990[ ). We used the brightness 



temperature of the 2_i — l_i E line in this calculation. Transitions between A and E species 
of CH3OH are strictly prohibited, therefore the total column density of the CH3OH molecule 
including both A and E species was estimated under the assumption that the abundance 
ratio between the A and E species is unity. The resulting column densities toward five YSOs 
are summarized in Table O 

The fractional abundance values of the CH3OH molecule at the protostellar positions 
were estimated by comparing the derived CH3OH column densities with the H2 column 
densities at the same positions. The H2 column densities were calculated using the 850 
/ im dust continuu m data observed with the JCMT SCUBA by the COMPLETE project 



( iRidge et al.l 120061 ). Under the assumption of optically thin dust emission, the H2 column 



density A(H2) is calculated by 

Cbeam 

N{^2) = ^ ^7;^, (13) 

where S"^*^^™ is the flux density per beam, f2beam is the beam solid angle, is the mass 
absorption coefficient per gram at 850 /im, and Bu(Tb) is the Planck function at a dust 
temperature Td. We assume the dust temperature of 12 K and kx = 0.012 cm^ g~^. The 
SCUBA image was smoothed to the 18" resolution in order to match the beam size of the 
CH3OH observations. The 850 /im flux densities, the derived column densities and abundance 
are summarized in Table [51 



The derived abundance ranges from 2.3 x 10"^ at SMM 1 to 9.4 x 10"^ at SMM 6, which 
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are in the same orde r of magnitude with t he CH3OH abundance in dark clouds, 2 — 5.1 x 10 ^ 
dPriberg et aD llQSsI : IPickens et aPboOoj ). The CH3OH abundance values toward SMM 6 
and SMM 11 are a few times larger than those toward the other SMM sources. These two 
sources with high CH3OH abundance are accompanied both with the CO outflow sho wn in 
the present observations and the shock-excited IR emission (jWalawender et al.ll2005bl ). On 
the other hand, there was no shocked IR emission found around SMM 2 and no significant 
outflow signature in our CO (1 — 0) data toward SMM 1 and SMM 3. Therefore, the higher 
CH3OH abundance values at SMM 6 and SMM 11 are considered to be influenced by the 
outflows. Such a variation of the CH3OH abundance is also found toward other YSOs. 



The CH3OH abundance was derived toward 30 Class and Class I sources (iBuckle fc Fuller 
2OO2I ). and the mean CH3OH abundance was 1.9 x 10"*^ for the sources with broad component 
whose FWHM is larger than 2.0 km s~^, and was 4.8 x 10~^ for the sources without broad 
component. The CH3O H abundance to^yard C lass sources is reported to be higher than that 
toward Class I sources (IBuckle fc Fullerll2002l ). however, our results show that the abundance 
depends on the presence of the broad line component rather than the evolutionary state of 
the sources. 

Next, we examine the CH3OH abundance in the outflow component around the position 
"a". At the position "a", the integrated intensity of the redshifted CH3OH wing emission 
(^LSR > 9.6 km s~^) is 1.1 K km s~^. Under the assumption of the LTE condition with 
an excitation temperature of 50 K, which is same as that of the CO emission, the column 
density of the CH3OH molecule was estimated to be 6.6 x 10^^ cm~^. Since dust continuum 
emission at the position "a" is too faint to estimate the H2 column density, we used the CO 
line intensity in the same velocity range as the CH3OH emission to calculate the H2 column 
density of the high velocity gas at this position. Under the assumption of the optically thin 
CO emission with an excitation temperature of 50 K, the CO column density was derived 
to be 1.8 X 10^^ cm-2. The abundance ratio of CH3OH to CO is [CH30 H/C0]= 3.6 x 10"^ 
and comparable to the abundance ratio measured in other outflows (I Gar ay et al.l |2002| ) . 
Assuming the CO abundance with respect to H2 is 10""^, the CH3OH abundance ratio, 
[CH3OH/H2], is 3.6 X 10^^. This value is more than 2 orders of magnitude larger than the 
abundance toward the protostars in Bl and the abundance in dark clouds. It is obvious that 
the CH3OH abundance is enhanced in the outflow. 



4.1.2. SiO 



The column density of the SiO molecule was estimated using Equation (1) under the LTE 
assumption. The rotation constant of the SiO molecule is 21,787 MHz (ILovas fc Krupenie 
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19741 ) and the dipole moment of the SiO molecule is 3.1 D. We assumed the optically thin 



condition and adopted an excitation temperature of 50 K, the same assumption used to 
calculate the CO and CH3OH column densities of the outflowing gas. Each SiO line proflle 
was divided into the "high-velocity" and "low velocity" components in order to examine the 
difference of the abundance relative to CO in these two velocity ranges. The threshold of the 
two ranges was set to be Vlsr = 9.6 km s~^ following the velocity criterion used to separate 
the CO line core and wing component. The derived SiO column densities at flve positions 
(SMM 2 Wl, W2, W3, E, and SMM 6) are summarized in Table El The column density 
of the SiO molecule was derived to be 1 — 2 x 10^^ cm~^ in the low-velocity component. 
The high velocity component was detected only at the positions of SMM 2 W2 and SMM 2 
W3. The SiO column density in the high-velocity component of SMM 2 W3 was 5.4 x 10^^ 
cm~^, which is a few times higher than the values measured in the low- velocity component. 
The relation between the total (high velocity + low velocity) SiO column density and the 
outflow velocity, which is determined from the terminal velocity of the SiO line, is shown in 
Fig. m The SiO column density tends to increase as the outflow velocity increases. Similar 
correlation between the SiO col umn density and th e outflow velocity was also found in other 
protostellar outflows studied by iGaray et al.l ( 12002| ). These results support the idea that the 
more SiO molecules evaporate from dust gr ains into gas phase by collisions of the grains with 
higher velocity molecules in the outflo ws (Schilke et al.l 19971 ). It should be noted that the 
SiO column density values reported in Garay et al.l ( 2002 ) are 1 order of magnitude lower 
than th ose derived here. Th is is probably because the region covered by the telescope beam 
used by iGaray et al.l ( 120021 ) with HPBW of 57" was twice larger than that covered by our 
38'.'4 beam. Since the size of the SiO emitting region is supposed to be smaller than the 
beams of the single-dish telescopes, the beam-averaged column density becomes lower if the 
source is observed with the larger beam. In addition, most of the sources in the sample of 



Garay et al.l (120021 ) are located further away in distance. 



The narrow (FWHM < 1.0 km s ^) SiO line as seen in SMM 6 was also detected in NGC 



1333 ( Lefloch et al.lll998l : ICodella et al 



19991 ). L1448-mm, and L1448 IRS3 (Ijimenez-Serra et al 



2004). The se authors proposed d ifferent mechanism for the SiO enhancement. iLefloch et al. 



( 119981 ) and ICodella et al.l ( Il999l ) suggested that the emission is arisen from the postshock 
equilibrium gas after the interaction o f prot ostellar jets with surrounding dense clumps. On 
the other hand, iJimenez-Serra et al.l ( 120041 ) proposed that the emission is produced by the 
shock precursors associated with the protostellar jets. The narrow-line component around 
SMM 6 was detected only toward S MM 6 and does no t have extended distribution like 
the one detected around NGC 1333 (ILefloch et al.lll998l ). The other difference is that the 
distribution of the narrow SiO line emission does not coincide with the high-velocity CO 
component in NGC 1333. These differences indicate that the mechanism to produce the 
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narrow SiO line emission is different in SMM 6 and NGC 1333 and it is likely that the shock 
precursor component produces the narrow line detected toward SMM 6. 

The abundance of SiO relative to CO was estimated by comparing the column densities 
of these molecules. Since the beam size of the SiO observations is larger than that of the 
CO, the CO map was smoothed in order to match the angular resolution with the SiO. 
However, our CO map is undersampled, and does not recover the entire flux. Therefore, 
we assumed that the spatial distribution of the CO flux between the observed points varies 
smoothly, and interpolated the CO flux values at the missing points using the neighboring 
data. This smoothing makes the derived CO column density ~ 25 - 40 % smaller than the 
value without smoothing. The CO spectra in the "low velocity" range contain the emission 
from the quiescent ambient material that is unr elated with the outflows . Therefore in the 
calculation we adopted a method suggested by iMargulis fc Ladal ( 119851 ) in which the line 
intensity originated from the low-velocity outflow is assumed to be the same as the line 
intensity at the boundary velocity between the line core and wing components; Vlsr = 2.1 
km s~^ and 9.6 km s~^ for the blue- and redshifted outflow components, respectively. The 
CO integrated intensities in the low velocity range are estimated by assuming the uniform 
intensity over the integrated velocity range, 4.0 — 9.6 km s~^ toward the positions along 
the SMM 2 outflow and 2.1 — 9.6 km s~^ toward SMM 6. The integrated velocity range is 
defined to cover all the SiO emissions at each position. The derived SiO abundance values 
are summarized in Table El Assuming [CO/H2] = lO^'*, the SiO abundance relative to H2 is 
10^'' — 10 ~^, 5-6 orders of magnitude higher than those measured in dark clouds (< a few 



xlO , IZiurvs et al.lll989l ). The highest SiO/CO abundance of 4.7 x 10 ^ was obtained in 



the high velocity component at SMM 2 W3. This abund ance is 3 orders of magnitude higher 
than the abundance estimated toward several outflows (IGaray et al.ll2002l ) and comparable 
to the abundanc e at the extremely h igh-velocity (EHV) jet component in the L1448-mm 
outflow (~ 10~^; iBachiller et al.lll99ll ). Such a high SiO abundance value and high terminal 
velocity (~ 50 km s~^, corrected for the inclination angle of 57?3) suggest that the high- 
velocity SiO emission observed at SMM 2 W3 also arises from the jet component. 



4.2. Outflows as an Engine of Interstellar Turbulence 



Since protostellar outflows are considered to be the most likely mechanism to maintain 
the turbulent energy of the cloud (ILi fc Nakamurall2006l ). we examine the turbulent energy 
budget of Bl main core following the method described in IStanke fc Williams! (120071 ). The 
energy-loss rate due to the turbulence decay, Lturb, and the energy input rate by the outflows. 



are compared. 
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According to numerical studies ( iMac Low fc Klesserul2004l ). the supersonic turbulence 
dissipates in a short period comparable to the free-fall timescale tg. Thus, Lturb is estimated 
by 



-'^turb 



(14) 



where EtuA = 3/2McorcC'"iD is the energy of the turbulence, ctid is the one-dimensional 
velocity dispersion. Since the observed region of the present s tudy is well matched with 
the extent of Bl main core identified with the NH3 observations (IBachiller et al.lllQQOl ). it is 
reasonable to adopt the mass M^ore and the one-dirnensio nal velocity dispersions (Tid of Bl 
main core (65 Mq and 0.43 km s~^; IBachiller et al.lllQQOl ) in this calculation. The free-fall 



time iff is calculated to be 3.5 x 10^ yr using a diameter of 0.5 pc and a mass of 65 Mq. 
Therefore, the turbulence energy loss rate is calculated to be Lturb = 8.5 x 1O^^L0. 

On the other hand, Lgam is calculated under the assumption that all the outflow mo- 
mentum is transferred to the turbulent motion of the cloud. 



1 



^gam 



(15) 



where Mc 



outflow 



/(JsD is the mass of the cloud which is accelerated up to the velocity 



C3D5 -^outflow is the total momentum input rate by outflows, and ctsd = y/ScriB is the three- 
dimensional velocity dispersion of the clump. The total momentum flux provided by the 
three outflows from SMM 2, SMM 6, and SMM 11 is Foutflow = 1.8 x IQ-^Mq km s-^jt-\ 
Therefore, the energy input rate is calculated to be Lgain = 1-1 x lO^^L©. This energy input 
rate is 1 order of magnitude smaller than the turbulence energy-loss rate. 



There are some uncertainties in estimating the outflow energy input rate. The first issue 
is the effect of optical depth of the CO emission. Since our calculation assumed the optically 
thin CO wing emission, the energy input rate derived here is considere d to be the lower 
limit. The mean opacity of the CO wing emission toward 16 YSOs is ~ 4 (jCabrit fc Bertout 



19921 ). and this value makes the energy input rate four times larger than that derived under 
the optically thin assumption. The second issue is the contribution of the low radial velocity 
component hidden in the line core. This effect becomes significant if the outflow axis is close 
to the plane of the sky, because most of the outflowing gas is moving along the transverse 
direction and shows the same radial velocity as the ambient gas. If the contribution of 
th e "hidden" low radial v elocity component is corrected by using the method described 
by iMargulis fc Ladal (Il985l ). which is also used to estimate the column density of the low- 
velocity outflow component in Section 4.1, the resulting energy input rate becomes twice as 
large as that calculated using the high- velocity wing component alone. The third issue is the 
contribution of the outflows that were not identified in the present observations. In addition 
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to the outflows from SMM 2, SMM 6, and SMM 11, otli er two outflows from S MM 1 and 
SMM 3 were identifled in the CO (3 — 2) observations by iHatchell et al.l ( l2007bl ). However, 
the contribution of these outflows to the total energy input is unhkely to be signiflcant; the 



momentum flux value of the SMM 3 outflow is only 1.8 x 10 ^M© km s ^yr ^ (IHatchell et al. 
2007bl ). which are more than a factor of 20 smaller than that of the SMM 2 outflow. In 



addition, the redshifted component of the SMM 1 outflow is treated as a part of the outflows 
from SMM 2 and SMMll in our calculations. In order to estimate the contribution of 
the weaker outflow components, we have integrated over all high-velocity component in the 
mapped area, and estimate the dynamical parameters of the high-velocity gas in this region. 
The total mass and momentum of the high velocity gas are derived to be 5.1 x 10"^ Mq and 
0.48 Mq km s~^, respectively. If we assume that the dynamical timescale is ~ 2 x 10^ years, 
which is a typical timescale of the outflows in this region, the energy input rate is calculated 
to be 1.4 X 10"'^ Lq. The energy input rate derived here is only a factor of 1.3 larger than 
that from SMM 2, SMM 6, and SMM 11 outflows, suggesting that most of the mechanical 
energy is supplied by the three major outflows. If we take the above issues into account, the 
outflow energy input rate can be 1 order of magnitude higher and becomes comparable to the 
turbulence energy decay rate. On the other hand, the efficiency of the energy transportation 
would be less than unity, because the extent of the outflow from SMM 2 and SMM 11 has 
already exceeded the dense part of the Bl main core region. Therefore, we conclude that 
the outflow energy input is not enough to maintain the turbulence in Bl main core. This 
result implies that the turbulence can not supply enough force against the gravity. 



The magnetic fleld in Bl was measured by lCrutcher et al.l (119941 ) to be 30 ± 4 yuG and 
16 ± 3 /iG in the NH3 core (r = lf5) and in the envelope (r = 9'), respectively, using OH 
Zeeman effect. The magnetic fleld strength, the size and the mean density of the NH3 core 



in Bl are well reproduced by the magnetically supercritical, collapsing core (ICrutcher et al. 



I994J ). This result indicates that the magnetic fleld is not able to support the core. Since 
the mass-to-flux ratio ( M/<l>) is indepen d ent o f the distance to the object, the difference of 
the distances assumed in ICrutcher et al.l (119941 ) and the present study would not essentially 
affect the conclusion. 

Since Bl is located in the Perseus cloud complex which contains several active star- 
forming regions, it is interesting to c o mpar e the star formation in Bl with that in other 
star- forming regions. |j0rgensen et al.l (120081 ) reported that there are eight protostars (i.e.. 
Class and Class I sources) out of nine YSOs in "Bl tight group", which corresponds to 
the Bl main core region. Since the fraction of protostars (8/9) in Bl is higher than the 
other star forming regions in Perseus, such as NGC 1333 (35/102) and IC 348 (11/121), they 
suggested that the star formation activity in Bl has recently initiated. In order to compare 
the turbulent energy budget in Bl with that of the more evolved star forming region, we 
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have calculated the turbulence supply rate and decay rate of the NGC 133 3 region. In the 
NGC 1333 region, six outflows were identifled by the CO observations by iKnee fc Sandell 
( 120001 ) ■ The total momentum injection rate from these outflows is ~ 5.6 x 10~^ Mq km s~^ 
yr~^, which is ~ 30 times larger than that in B l. Using the velocity dispersion of 0.85 km 
s~^, which is derived from the NH3 observations (jSchwartz et al.lll978l ). the energy input rate 
Lgain is estimated to be 6.8 x 10~^ Lq in NGC 1333. On the other hand, if we adopt the mass 
of the core to be 319 M©, the turbulent decay rate Ltmh is estimated to be 1.4 x 10^^ Lq. 
These two values agree within factor of 2, thus the out fl ows a re powerful enough to maintain 
the turbulence in NGC 1333. Since iKnee &: Sandelll (I2OOOI ) assumed the lower excitation 
temperature of 20 K and did not apply the inclination correction, the energy input rate is 
underestimated compared to the case using our assumption of the excitation temperature 
of 50 K and the inclination angle of 57?3, in which the outflow mornentum flux becomes a 
factor of 6 larger than the one originally derived by lKnee fc Sandelll (120001 ). This makes the 
outflow energy input superior to the turbulence energy loss. 



The magnetohydrodynamical simulation (INakamura fc Lil 120071 ) shows that the cloud 
contraction proceeds in the initial phase of the cluster formation due to the decay of turbu- 
lence and decrease of the cloud momentum. This is consistent with the observed picture of 
Bl main core. The cloud reaches an equilibrium state in the later phase of the simulation 
because of the energy input by the outflows. As their simulations, when the outflows from 
the young Class sources in Bl have enough evolved, it might be possible that the outflows 
have a considerable contribution to the energy balances of the turbulence. 



5. SUMMARY 

In order to study how outflows from protostars influence the physical and chemical 
conditions of the parent molecular cloud, we have observed Bl main core, which harbors 
four Class and three Class I sources, in the CO (J = 1 — 0), CH3OH {Jk = 2k — Ik), 
and the SiO {J = 1 — 0) lines using the NRO 45 m telescope. We have identifled three 
CO outflows in this region; one is an elongated (~ 0.3 pc) bipolar outflow from a Class 
protostar Bl-c in the submillimeter clump SMM 2, another is a rather compact (~ 0.1 pc) 
outflow from a Class I protostar Bl IRS in the clump SMM 6, and the other is the extended 
outflow from SMM 11. No signiflcant outflows from other YSOs were identified in our CO 
( J = 1 — 0) map. The CH3OH emission lines are distributed around Bl main core with two 
peaks; one peak is located at the position of SMM 6, toward which the CH3OH emission 
in the line core velocity (5.5 < Vlsr < 7.5 km s~^) is enhanced, and the other is 120" west 
of SMM 2, at the head of the redshifted CO lobe of the SMM 2 outflow, toward which the 
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CH3OH line shows broad redshifted wing emission with a terminal velocity of 13 km s~^. 
In the redshifted lobe of the SMM 2 outflow, the SiO line also shows broad redshifted wing 
whose terminal velocity reaches up to 25 km s~^. 

The CH3OH and SiO abundances are significantly enhanced at the redshifted lobe of the 
SMM 2 outflow. It is likely that the shocks caused by the interaction between the outflow 
and ambient gas enhance the abundance of SiO and CH3OH in the gas phase. Toward 
SMM 6, the CH3OH and SiO lines are narrow and peak at the systemic velocity, however, 
the abundances of these molecules are enhanced compared with those in the quiescent dark 
clouds. Since the SiO emission around SMM 6 is localized at the position of the protostar, 
the origin of this narrow component is likely to be the shock precursor rather than the fossil 
outflow. 

The turbulence energy decay rate and the outflow energy input rate are compared to 
investigate the energy budget in Bl main core. It is found that the energy input rate, 
1.1 X 1O~^L0, is approximately 1 order of magnitude smaller than the turbulence decay 
rate of 8.5 x 10~^Lq. Even though several uncertainties that increase the energy input rate 
are taken into account, the outflows in Bl main core are not energetic enough to sustain 
the turbulence in this region. Since t he magnetic fleld is ri ot strong enough to prevent the 



gravitational collapse of Bl main core (ICrutcher et al.lll994l ). it is likely that Bl main core is 



gravitationally unstable. The high fraction of protostars among YSOs (8/9) in Bl main core 
suggests that the star formation activity in this region has been started very recently. In the 
more evolved star-forming region, NGC 1333, the cloud can be supported by the turbulence 
driven by th e outflows. These resu lts are consistent with the picture of the cloud evolution 



proposed by iNakamura fc Lil (120071 ) , in which molecular clouds contract due to the decay of 
turbulence in the early evolutionary phase, and become stable once the energy input from 
protostellar outflows become sufficient enough to balance the energy loss. 
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Table 1. Young Stellar Objects in Bl 



Name 


a 


5 


-^bol 




Class 


other names 




(J2000) 


(J2000) 




(K) 






SMM 1 


03 33 21.5 


31 07 29 


<2.5 


<25 





Bl-bN and Bl-bS 


SMM 2 


03 33 18.0 


31 09 32 


3.7 


53 





Bl-c 


SMM 3 


03 33 16.3 


31 06 53 


2.6 


32 





Bl-d 


SMM 6 


03 33 16.6 


31 07 47 


1.3 


158 


I 


Bl-a, Bl IRS, IRAS 03301+3057 


SMM 11 


03 33 27.3 


31 07 08 


0.4 


117 


I 




LkHa 327 


03 33 30.4 


31 10 51 






I 


IRAS 03304+3100 



Note. — Positions of the SMM sources are referred from lWalawender et al.l ( l2005bl ). and that 
of L kHa 327 is obtained fr om Simbad. Lboi, Tboi, and Clas s of the SMM sour c es are referred 
from iHatchell et al.l ( l2007al ). The class of LkHa 327 is from lWalawender et al.l ( l2005bl ). 



Table 2. Observed Molecular Lines 



Molecule Transition Frequency (GHz) (K) 



CO 


J = 


-- 1 - 







115.2712 


5.53 


CH30H 


Jk 


= 2_ 


1 - 


- l_i E 


96.73939 


4.64 




Jk 


= 2o 




lo 


96.74142 


6.79 




Jk 


= 2o 






96.74458 


12.2 


SiO 


J = 


= 1 - 







43.42386 


2.08 



Note. — The upper state energies of the CH3OH E and A'^ 
symmetric state are the values relative to the ground state of 
l_i E and Oq A'^, respectively. The l_i E level is 7.9 K above 
the Oo A+ level. 
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Table 3. CO Line Intensity Ratio 



Source 


Position (CO 3-2) 


Integ. Vel. Range 


/ TMB,3-2dv 


/ TMB,i-odv 


T 




a(J2000) (5(J2000) 


(km s~^) 


(K km s-^) 


(K km s-^) 


(K) 


SMM 2 red 


03 33 15.866 31 09 49.13 


10 - 20 


24.8 


4.5 


46 


SMM 6 


03 33 16.257 31 07 49.13 


-10 - 2 


33.2 


6.2 


42 



Table 4. Outflow Parameters 



Component 


M 




P 




M 






F 






(Me) 




(Mgkm s 




(Me yr-i) 


(yr) 


(Mekm^ s-2) 


(Mekm s"^ 




SMM 2 


Blue 


1.5x10" 


-2 


0.14 




4.8x10"^ 


2.9x10^ 


0.66 


4.8x10" 


-6 


Red West 


5.3x10" 


-3 


3.9x10" 


-2 


1.8x10'^ 


2.1x10^ 


0.17 


1.8x10" 


-6 


Red East 


5.2x10" 


-3 


3.3x10- 


-2 


2.1x10-^ 


1.6x10^ 


0.11 


2.1x10" 


-6 


Total 


2.5x10" 


-2 


0.21 




8.7x10-^ 




0.93 


8.7x10" 


-6 


SMM 6 


Blue 


1.2x10" 


-3 


1.2x10" 


-2 


1.1x10"^ 


1.1x10^ 


5.6x10-2 


1.1x10" 


-6 


Red 


2.3x10" 


-3 


1.7x10" 


-2 


9.8x10"^ 


1.8x10^ 


7.0x10-2 


9.8x10" 


-7 


Total 


3.5x10" 


-3 


2.9x10" 


-2 


2.0x10"^ 




0.13 


2.0x10" 


-6 


SMM 11 


Blue 


1.4x10" 


'2 


0.14 




6.1x10"-^ 


2.3x10^ 


0.72 


6.1x10" 


~6 


Red 


5.1x10" 


-3 


3.0x10" 


-2 


9.4x10-9 


3.2x10^ 


8.7x10-2 


9.4x10" 


-7 


Total 


1.9x10" 


-2 


0.17 




7.0x10-^ 




0.81 


7.0x10" 


-6 
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Fig. 1.— CO (J = 1 - 0) line profiles toward (a) 120" west of SMM 2, (b) SMM 6, (c) SMM 
1 and (d) 5^5 southwest of SMM 6. The velocity resolution is 0.3 km s~^. Vertical lines 
indicate the outflow velocity boundaries. 
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3^3;4Q S3'3Q 33:2Q 33:10 33:00 

Right Ascejisjon (J2D00} 

Fig. 2. — Distribution of the high velocity component of the CO (J = 1 — 0) hne. The 
blueshifted (—8.1 km < Vlsr < +2.1 km s~^) and the redshifted (9.6 km s^^ < Vlsr < 
19.8 km s~^) components are drawn with blue and red contours, respectively. The first 
contour level and the contour interval are 3cr, which correspond to 0.5 K km s~^ for blue and 
0.68 K km s~^ for red contours, respectively. The gray-scale image in the background is the 
Spitzer IRAC2 4.5 /xm image. Circle an d cross marks indicate the positions of protostars 
and H2 knots ( IWalawender et al.ll2005bl ). respectively. X marks indicate the positions "a" 
to "d" in FigJD 
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Fig. 3.— CH3OH {.Jk = 2o - lo = 2o - lo A, and 2_i - l_i line profiles toward (a) 
120" west of SMM 2, (b) SMM 6, and (c) SMM 1. The velocity resolution is 0.3 km s'^. The 
reference frequency in the velocity calculation is set to be 96.73939 GHz, the rest frequency 
of 2_i — l_i E transition. 
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Fig. 4. — Integrated intensity map of the CH3OH lines. Black contours represent the inte- 
grated intensity of the Jk = 2_i — l^i E and 2o — lo A'^ lines {A + E map). Red contours 
show the integrated intensity of the 2o — Iq E line. Both the lowest contour level and the 
interval are 0.23 K km s~^ {6a) for the A + E map. The lowest contour level is 0.13 K km 
s~^ and the contour interval is 0.043 K km s~^ (la noise level) for the 2q — 1q E line. The 
integrated velocity range is —3.6 < Vlsr < 20.4 km s~^ for the A + E map and 2.9 < Vlsr < 
10.2 km s~^ for the 2o — lo line. The position "a" and the protostars are marked with a 
cross and circles. 



- 30 - 



I I I ■■ t I ■ I I I I I ■ I I I ■ I I I I I ■ I I I I I I I I I ■ I q ■ I I I I I I ■ I I 




3:33:40 33:30 33:20 33:10 33:00 

Right Ascension (J2000) 



Fig. 5. — Integrated intensity map of the CH3OH {Jk = 2_i — l_i E) line in the hne core 
velocity range (5.5 km s^^ < Vlsr < 7.5 km s~^, black contours) and the higher velocity 
range (2.5 km s~^ < I^lsr < 5.5 km in blue contours and 7.5 km s~^ < Vlsr < 20 km s^^ 
in red contours). Note that the blueshifted component around the position "a" (black cross) 
is due to a contamination of redshifted emission of the Jk = '^o ~ lo line. The lowest 
contour and contour intervals are 3a noise level of each map, 0.12 K km s~^ for the line core, 
0.09 K km s~^ for the blueshifted, and 0.24 K km s~^ for the redshifted component. 
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Table 5. CH3OH Column Density and Abundance 



Source 


/ TuBdv 

(K km s-i) 


N 

(cm-2) 


(Jy l")eam^"^) 


A^(H2) 

(cm-2) 


[CH3OH/H2] 


Shocked IR 


SMM 1 


1.9 


2.3x10^^ 


1.0 


1.0x10^3 


2.3x10"'' 


No 


SMM 2 


1.7 


2.0x10^^ 


0.87 


8.6x10^2 


2.4x10-*^ 


No 


SMM 3 


1.3 


1.6x10^^ 


0.64 


6.3x10^2 


2.5x10-'' 


Yes 


SMM 6 


3.4 


4.1x10^4 


0.44 


4.4x10^2 


9.4x10-'' 


Yes 


SMM 11 


0.58 


7.1x10^=^ 


0.12 


1.1x10^2 


6.2x10-'' 


Yes 



m» ji|H4 . 



3:;J3:40 33' 30 33;£0 33:10 33iD0 

Rj^t Ascensian (J2(K>0} 



.^i-ir.«r 



Fig. 6. — SiO {J — 1 — 0) observing positions overlaid on the Spitzer IRAC2 image. The large 
circles with solid and broken lines indicate the positions where the SiO line was detected and 
not detected, respectively. The size of the circles corresponds to the HPBW in the frequency 
of the line (43.4 GHz). 



Table 6. SiO Column Density Along the Outflows 



High Velocity Low Velocity 

« 5 jTMB,siodv 7V(SiO) [SiO/CO] jTMB,siodv N{SiO) [SiO/CO] 

Source (J2000) (J2000) (K km s-^) (cm'^) (K km s-^) (cm-^) 



SMM 2 W3 


3 


33 


09.3 


31 


09 57.9 


2.6 


5.4 X 




4.7 X 10"^ 


0.78 


1.6 


X 




9.1 


X 


10- 


-4 


SMM 2 W2 


3 


33 


12.2 


31 


09 49.2 


0.69 


1.4 X 




1.4 X 10"^ 


0.83 


1.7 


X 




1.0 


X 


10" 


-3 


SMM 2 Wl 


3 


33 


15.1 


31 


09 40.6 










0.47 


9.8 


X 


10^2 


5.1 


X 


10- 


-4 


SMM 2 E 


3 


33 


29.7 


31 


08 57.4 










0.59 


1.2 


X 


10^3 


2.7 


X 


10- 


-4 


SMM 6 


3 


33 


16.6 


31 


07 46.8 










1.1 


2.3 


X 


10^3 


2.1 


X 


10- 


-3 
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Fig. 7. — SiO {J — 1 — 0) line profiles along the outfiow axis of SMM 2 and at the position 
of SMM 6. The velocity resolution is 0.5 km s~^. 
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oJHlDiir valodfy (km sr^) 

Fig. 8. — Relation between the outflow terminal velocity measured from the systemic velocity 
at each point and the column density of the SiO molecule. The errors are calculated from 
the absolute intensity calibration uncertainty (~ 15%). 



